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multiplied by uo(y) = e * y2 at infinity. Consequently, we need the series to terminate, 
which requires e m = 2m + 1 for some m. Thus, 



^ *. ( 1 
E n = —e n = huj\n+- 



HO energy levels (18-1) 



Quantized energy levels of a harmonic oscillator. The ground state (zero- 
point) energy is E = ^hu, the energy levels are equidistant. 

Note. This feature allows us to identify the HO not only with a particle in potential 
V(x) = ^mu 2 x 2 , but also with a system of noninteracting (bosonic) particles. 

Therefore, a mode of an electromagnetic field of frequency u can be viewed as a HO 
with frequency u; n photons in that mode correspond to the n-th occupied state 
of the HO. The uncertainty in x and p of the HO ground state corresponds to the 
"vacuum fluctuations" of the electromagnetic field (x) = 0, (re 2 ) ^ corresponds to 
(E) = 0, (E 2 ) ^ etc. 

For given e = 2n + 1, the recursion relation (*) 

(m + l)(m + 2)c m+2 = (2m - e n + l)c m = (2m - 2n)c m (18-2) 

yields 

• c 2 = -f^c 

. Ca - ^2n ro _ (4-2n)(0-2n) 

• C 4 — 3 4 C 2 - ^2.3.4 C 

and in general 

. n , k n(n-2)---(n-2k + A)(n-2k + 2) 
c 2k = (-2)^ (2fc)| ^c , 



< 2k < n, n even (18-3) 



for the even coefficients. 

For the odd coefficients we have 

• C3 = ^ Cl = (-2)^ Cl 

_ 6-2n _ / n\n-3 _ ( q\2 (n-3)(n-l) 

• C 5 - ^5" C 3 - l-^J^- c 3 - l-^J 1.2-3-4-5 Cl 

and in general 



C 2 fc+1 



/ , fc (w-l)(7i-3)-(n-2fc + 3)(n-2fc + l) CT 



(2fc + l)! 



< 2k + 1 < n, n odd. (18-4) 
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by 



The eigenfunction u n (x) for energy level n with energy E n = hcu (n+ |) is given 



n/2 



n=0 
(n-l)/2 



U 



n (y) = e h/y ^2 C ^+lV 



2fc+l 



for even n, 



for odd n, 



(18-5) 



(18-6) 



n=0 



The coefficient Co or c\ has to be chosen such that the wavefunction is normalized, 
and y is related to the position coordinate x via y = x /^x. The quantity \ — 

a r a V h ^ J y 

has units of length and defines the natural quantum length scale for the harmonic 
oscillator. Apart from the normalization, the polynomials 



n/2 

My) = $^ c 2fc2/ 2fe 

fc=0 

(n+l)/2 

Ms/) = c ^+iy 2k+l 

k=0 



is-?) 



(18-8) 



are the Hermite polynomials H n (y). The Hermite polynomials obey the following 
relations: 



H n "(y) - 2yH n \y) + 2n# n (y) = 
H n+1 (y) - 2yH n (y) + 2nH n ^(y) = 
H n+1 (y)-H n '(y) + 2yH n (y) = 



(defining equation) 



X>n(!/) 



n=0 



!T B (y) = (-l)V 



alternative definition \ 
of Hermite polynomials J 

i 2 i 



(18-9) 
(18-10) 
(18-11) 

(18-12) 
(18-13) 



Since the wavefunction belonging to level n is u n (y) = Ce | H n (y), in order to 
normalize it, we need to calculate 

/OO PCX) 
dy\u n (y)\ 2 = / dy\C\ 2 e~y 2 H n \y) (18-14) 
-OO «/ — OO 

The Hermite polynomials are real. One can show that 



/OO 
dye-y 2 H n \y) = 2 n n\^ 
-OO 



18-15) 
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Normalization actually requires dx\u n (x)\ 2 = 1, but since x and y = ^J^x are 
related by a constant factor, dx\u n (x)\ 2 — J ^ dy\u n (y)\ 2 . 



(18-16) 
(18-17) 
(18-18) 
(18-19) 
(18-20) 
(18-21) 

Consequently, the lowest eigenfunctions look like (Fig. I). The eigenfunctions of the 



H (y) 


= 1 


Hi(y) 


= 2y 


H 2 (y) 


= V - 2 


H 3 (y) 


= % 3 - 12y 


H*(y) 


= 16y 4 - A8y 2 + 12 


H 5 (y) 


= 32y 5 - 160y 3 + 1202/ 



\ A 




U ■— -"V V— 




\— : 
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> 

Figure I: HO eigenfunctions. 



HO look the same in momentum space, since the Hamiltonian is symmetric in x and 
p, and 



x = ih 



d_ 
dp 



,_hd 
^ i dx 



p = p 



in position space ip(x) 
in momentum space <j)(p) 



(18-22) 
(18-23) 
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The harmonic oscillator ground state, being a Gaussian function with no spatial de- 
pendence of the complex phase, has minimum uncertainty allowed by the Heisenberg 
relation: 



AxAp > 



h 

2 
Ti 



for the ground state 
for any excited state 



(18-24) 
(18-25) 



• Show Bose-Einstein condensate expansion 

• Thermal cloud 

• isotropic expansion for anistropic trap 



• v 

■ 2m 



\kT) 



condensate: 

pi (x J_ ^ m 

2m xi) 2 h 



anisotropic expansion 



HO: operator method 

There is an elegant and instructive way to derive the HO eigenstates without directly 
solving the SE. Instead, we use commutation relation between operators. We start 
by writing the Hamiltonian in dimensionless form 



H = hu 



= Tiui 



p 



+ 



moo 



Imhuj 2H 



p 



-x 



\j2mhuj 

2 



X 



+ 



x 

Xq 



(18-26) 
(18-27) 

(18-28) 



with p 2 = 2mhuj, Xq 1 



2h 

mui ' 



Classically we can write, 



H d = huj 



x . p \ I X 
,^0 Po) \x Po 
however, since in QM p and x do not commute, we have 



P 



\xq Po 



x p 
x po 



X 

x 

= (— 



+ 



x po 



[x,p] . 



(18-29) 

(18-30) 
(18-31) 
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/raw 1 



Using the commutator [x,p] = thgp - pihf p = ih and ^ - y 2h - 2h 

have 



we 



H — ftuj 

= Tiu 
= %ui 



x 
x 



x . p \ I x . p 



xo Po/ \%o Po 
x . p\ ( x . p 



x p 



[x,p] 



— I- 



+ 1- 



Po/ \xo Po 
We can define a new, non-Hermitian operator by 

x . p 

a := h i — 

x po 

Consequently, the Hermitian conjugate operator is 



(18-32) 

(18-33) 
(18-34) 



(18-35) 



x p 
xo Po 



t 



x . p 
x p 



(18-36) 



since ft — p, x^ — x 



Note. The Hermitian conjugate operator 0t of any operator is defined by the relation 



dxi> 2 \x)0^ l {x)= / dx (Oi> 2 (x))* ^(x) 



(18-37) 

for any well-behaved wavefunctions ipi(x), ip 2 (x). 

Consequently, for any operator O = C\0\ + c 2 2 , where C\ and c 2 are complex 
numbers, 

6 f = (ciOi + c 2 6 2 ) ] = c^Ol + c 2 *6\ (18-38) 
and for any operator O = 0\0\ we have Cfi = (0\d 2 ^j = 0\d\ 
Proof. See problem set. 



□ 



Using the operators a,aJ , we can write the Hamiltonian for the HO in the partic- 
ularly simple form 



H = %uj ( aJa + 



1 



48-39) 
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Rather than being explicitly defined in terms of x, p, . . . , and operator can be defined 
through its commutation relations with other operators. Let us look at a, a): 



[d,d ] ] = 



x 
x 



p X 

I—, — 

Po X 



■ P 

I — 

Po 



= I 



p X 


— i 


x p 


Po' Xo_ 




x ' p 0m 



Pox 
% 

i h 
h i 

1, 



[x,P]) 



(18-40) 
(18-41) 
(18-42) 

2[p,x] (18-43) 

(18-44) 
(18-45) 

where we have used [x, x] = = [p, p] . So we have 

[a,a f ] = 1 (18-46) 

[a,a] = [tf,a!]=Q (18-47) 

As will be elaborated on in 8.06, this defines a commutation relation for bosonic 
(quasi)-particles, i.e. particles whose wavefunction is symmetric under the exchange 
of two particles. For the commutators with the Hamiltonian, we have 



H, a J = hw [a^a, a] 



a)aa 



da) a] 



and 



H,a) 



— hu ( 

= hw (a'aa — (l + a) a) a) 
= —hud 



(18-48) 

(18-49) 
(18-50) 
(18-51) 



hid (a) da) — a) a) a) 



fiuj (a^acJ — (acJ — l)) 
hujd) 



(18-52) 

(18-53) 
(18-54) 
(18-55) 

We are now in the situation to calculate the spectrum of eigenenergies of the HO 
simply using those commutation relations. Let us first note that since H is quadratic 
in x and p, all eigenvalues must be positive: 

(E) = (A) = (T) + (V) (18-56) 



1 

2m 



dpcjf (p)p 2 <fi{p) + -muo 2 I dxil)*(x)x 2 il)(x) > 



18-57) 
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for any wavefunction ip(x) and its Fourier transform <f)(p). Before determining the 
eigenspectrum, let us define a convenient notation. 



State vector notation (Dirac notation) 

We have already argued that a physical state, (i.e., a physical system whose initial 
conditions have been prepared to the maximum extent allowed by QM), is described 
by a vector in an abstract vector space (Hilbert space), and that a wavefunction 
in position space is only one possible representation of the state. Alternatively, the 
state can be described in the momentum representation (wavefunction in momentum 
space), or by specifying the expansion coefficients when expanding the basis of energy 
eigenstates. Using a notation introduced by Paul Dirac, one of the creators of QM, 
we write the state as 

(18-58) 

and define 

: = J dx<p*(x)ip(x) (18-59) 
for any two states |0) whose wavefunctions are by ip(x), (j>(x). Dirac introduced 



bra -c- ket 



(18-60) 



So is called a "ket", and (0| a "bra". You can think of the "bra" as the transpose 
of the "ket" vector 



(yi v 2 • • •) 



c-number (complex number), 



18-61) 



but possibly for infinite-dimensional vectors. In this sandwich or Dirac notation, the 
expectation value of any operator A is given by 



Aj = (ip\A\ijj^ = / dxip*(x)Aip(x 



In Dirac notation, 

(m = 



dxip*(x)(f)(x) 



dx(fr* (x)ip(x) 



my 



An operator A acting on a state produces another state, symbolically 



18-62) 



(18-63) 



(18-64) 
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Consequently, 

(<P\A\^ = = j dx(j)*(x) (M{x)} = J dx<p*(x)Ai)(x) (18-65) 

The Hermitian conjugate operator A' is defined by 

(<j>\A^ = J dx (A(j)(x))* ^j(x) = (^j dx^j*(x)(A(j)(x))^ (18-66) 

= ($\A<i>)* (18-67) 
= (A(f)\ip) (18-68) 

In Dirac notation, the orthonormality condition for eigenstates \n), \m) reads 

(n\m) = J dxu n * (x)u m (x) = 5 nm , (18-69) 

where the expansion coefficients are 

Cn = J dxu n *i){x) = (n\ip) . (18-70) 

A bracket like (a\b) is a complex number, but a ketbra like \b) (a\ is an operator since 
acting on a state it produces another state 

• (18-71) 

state c-number 

One can show that the sum over all eigenstates J2 n \n) {n\ oi & Hermitian operator is 
the unity operator 

" "~ n ) (n\ = 1, (18-72) 



E 

n 

and 



,,) = i m = E ,„> m = E c„ ,„> ( — ) , 18 -73) 

Back to the operator treatment of the HO: Let us assume that we have found an 
energy eigenstate with eigenenergy E and let us denote that state by \E). Let us 
define a new state \ip) by having the operator a act on \E), \ip) := a\E). What 
happens if we act with the Hamiltonian 

H\iP) = Ha\E} (18-74) 

= (H,a + a#) \E) (18-75) 

= (-hcua + aE) \E) (18-76) 
= (E-hu)a\E) (18-77) 
= (E — huj) (18-78) 
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H,a 



—hua for the commu- 



Here we have used the previously calculated result 

tator, and the fact that complex (here real) numbers commute with everything. The 
above formula signifies that \ip) is also an energy eigenstate, but with lower energy 
E = —huj. Since starting from any eigenstate \E) we can repeat the procedure any 
number of times, 

a n \E) = \E- nhu) , (18-79) 



and the eigenenergy has to remain positive, (we have shown {^j}\H\ipj > for any 
state), there must exist a state |0) such that 

a|0) = (18-80) 

i.e., a state whose energy cannot be lowered further. 

Note. It is important to distinguish between |0) (lowest energy eigenstate, vector in 
Hilbert space) and (zero of the Hilbert space, vector of zero length). 

Note. Nothing implies that the state |0) has zero energy, in fact, 

H\0) = hu)(tfa+\)\0) = lhuj\0), (18-81) 
y 2 J a|o)=o 2 

so the ground state has eigenenergy Eq = ^hcu, this is the zero-point energy. In the 
context of identifying a HO at frequency u with an electromagnetic mode at frequency 
u), the ground state |0) is also called the vacuum (ground state has no excitations, 
photon number is zero): the vacuum has finite vacuum energy E = \Tiujq. 

What happens if acts on ground state? Let us define 

|l) := a 1 " |0> . (18-82) 

The tilde ~ is there to remind us that this state is not necessarily normalized, even 
if |0) is chosen to be normalized. 
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